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Abstract. The Kac-Ward formula allows to compute the Ising parti- 
tion function on any finite graph G from the determinant of 2 2fl ma- 
trices, where g is the genus of a surface in which G embeds. We show 
that in the case of isoradially embedded graphs with critical weights, 
these determinants have quite remarkable properties. First of all, they 
satisfy some generalized Kramers-Wannier duality: there is an explicit 
i— <*' equality relating the determinants associated to a graph and to its dual 

graph. Also, they are proportional to the determinants of the discrete 
critical Laplacians on the graph G, exactly when the genus g is zero or 
one. Finally, they share several formal properties with the Ray-Singer 
9-torsions of the Riemann surface in which G embeds. 

1. Introduction 

Most of the exact results for the two-dimensional Ising model rely on the 
so-called Pfaffian method. The idea, due independently to Hurst-Green [23] . 
Kasteleyn [26] and Fisher [18l [19] , is to associate to the given graph G an 
auxiliary graph Tq such that the dimer partition function on Tq is equal 
to the Ising partition function on G. The dimer model technology can 
then be applied to solve the Ising model on G. In particular, if G can be 
embedded in an orientable surface of genus g, then so can Tq, and its dimer 
partition function can be computed as an alternated sum of the Pfaffians 
of 2 2fl well-chosen skew-adjacency matrices [221 EEl E3 HO]. Even though 
this method is very standard at least for planar graphs (see for example 
the classical book [34] ) . it is not a very natural one. First of all, this so- 
called Fisher correspondence G i— > Tq is by no means unique: for example, 
each of the articles [26j [TOJ [TOJ El [Ml IE] contains a different version of it. 
Secondly, whatever the chosen correspondence, virtually all the geometric 
and combinatorial properties of G will get lost when passing to Tq - the 
only obvious exception being the genus. This will not be a problem if one is 
studying a topological class of graphs, such as all planar graphs, or all finite 
graphs of a given genus. However, if one is interested in a geometric class - 
as we will be - the Pfaffian method leads to unnecessary complications. 

There is another combinatorial method to solve the two-dimensional Ising 
model which, although naturally related to the Pfaffian method (see [SI 
Subsection 4.3]), is in our opinion much more natural. It is due to Kac- 
Ward [2H] (even though a rigorous proof awaited many years [TO]), was 
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originally formulated for planar graphs, and recently extended to any finite 
graph [32l IB]- With this method, no auxiliary graph is needed: the Ising 
partition function on a finite graph G is computed as an alternated sum of 
the square roots of the determinants of 2 2g Kac-Ward matrices naturally 
associated to the graph G. 

In this paper, we initiate the study of the critical Ising model on graphs 
of arbitrary genus, and we do so using these Kac-Ward matrices. What we 
mean by "critical Ising model" will be formally defined and thoroughly mo- 
tivated in Section [21 but for now, we will content ourselves with an informal 
description. Consider a finite number of planar rhombi {R e } e £E of equal 
side length, each rhombus R e having a fixed diagonal e and corresponding 
half-rhombus angle 6 e S (0, n/2), as illustrated below. 




Paste these rhombi together along their sides so that extremities of diag- 
onals are glued to extremities of diagonals. The result is a graph G with 
edge set E(G) = E embedded in a flat surface E with so-called cone-type 
singularities in the vertex set V(G) of G, and in the middle of the faces of 
G C E. We shall say that G is isoradially embedded in the flat surface E 
(see Definition [T]) . Using the high-temperature representation, we define the 
partition function for the critical Ising model on such a graph G by 

7e£(G) e£7 

where £(G) denotes the set of even subgraphs of G, and the critical weights 
v e are given by v e = tan(0 e /2). This is a natural generalization of the critical 
Z-invariant Ising model [I], which corresponds to the special case where E is 
a domain in the (flat) plane. Note that only a specific class of planar graphs 
admit a Z-invariant Ising model (see |291lll|). On the other hand, any finite 
graph can be isoradially embedded in a flat surface as explained above, and 
therefore admits a critical Ising model (Proposition 12. 1 | ) . 

Since its introduction by Baxter, the critical Z-invariant Ising model has 
been extensively studied, as well as its analog on the flat torus (see for 
example the papers [6l [7], where Boutillier and de Tiliere make use of the 
Pfaffian method). On the other hand, very little is known about the critical 
Ising model on graphs in higher genus. It is widely believed that such models 
should be discrete analogs of some conformal field theory pQ, but in genus 
g > 2, such statements are only supported by numerical experiments on 
very specific examples [12], [13] . It is our belief that, in order to try to tackle 
such outstanding conjectures, the Kac-Ward method will prove useful. 

Before doing so, we need to settle some fundamental questions about the 
critical Ising model on graphs of arbitrary genus and the associated Kac- 
Ward matrices, and this is exactly what the present paper is about. 
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In the general case of an arbitrary graph embedded in a topological sur- 
face, the (generalized) Kac-Ward matrices can be quite complicated (see 
Definition 1 in [8]). The first nice surprise is that for a weighted graph 
(G, x) embedded in a flat surface, the corresponding Kac-Ward matrices 
take a remarkably simple form. For each homomorphism <p from the fun- 
damental group 7Ti(£) of E to the group S 1 , we get a ip -twisted Kac-Ward 
matrix of order 2\E(G)\ whose determinant we denote by t^G, x) (see Def- 
inition [3]). For a specific type of 93's - the ones that belong to the set S of 
discrete spin structures on G C E - t v (G, x) turns out to be the square of 
a polynomial in the variables {x e } e . The main theorem of [8] then easily 
implies the following result (see Theorem 13.31 for the complete statement). 

Theorem. If all cone angles are odd multiple of 2ir, then the Ising partition 
function on the weighted graph (G, x) is given by 



Aes 

where g is the genus o/E and Arf(A) E Z2 the Arf invariant of the discrete 
spin structure A. 

Let us now turn to our main results. In a nutshell, we show that when the 
weight system is critical the corresponding Kac-Ward determinants t </3 (G, v) 
exhibit several remarkable properties. 

Firstly, these determinants turn out to admit a relatively simple combi- 
natorial interpretation, as described in Proposition 14.11 Furthermore, they 
satisfy the following duality property. 

Theorem. Let G be a graph isoradially embedded in a flat surface E, and let 
v be the critical weight system on G. The dual graph G* is also isoradially 
embedded in E, and therefore admits a critical weight system u*. If all cone 
angles are odd multiples of 2ir, then for any if, 



2 lV(G * )l II (l+cos(M)^(G">*) = 2l^ G )l J] (l+cos(0 e )) r^(G». 



This can be interpreted as a generalization of the celebrated Kramers- 
Wannier duality [31] from the case of planar graphs to the case of graphs of 
arbitrary genus. 

Our final result relates the Kac-Ward matrices to some a priori totally 
different operator. Given a weighted graph (G, x) and a homomorphism 
ip: vri(G) — > S , the associated discrete Laplacian on G is the operator 
A v = Af(G, x) acting on / E C v ^ by 



the sum being over all oriented edges e of the form (v, w). Note that if G is 
a (planar) isoradial graph, then the corresponding critical weights are given 
by c e = tan(# e ) (see [2Z])- We prove: 

Theorem. Let G C E be an isoradially embedded graph in a flat surface, 
and let us assume that all cone angles "& v of singularities v E V(G) are odd 




e*&E(G*) 



e£E(G) 



(A«7)(i>)= E Xe(f(v)-f(wMe)) 
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multiples of 2ir. If the genus of S is or 1, £/ien for any (p: 7Ti(S) — > S , 
t*(G,v) = (-1) n 2-xW J] '° S( ^, detA^(G,c), 

eeE(G) 1 + C ° S(0e) 

where N is the number of vertices v G ^(G) suc/i i/iai '& v /2tt is congruent to 
3 modulo 4, and x(^) = |y(G)| — |J5(G)|. On the other hand, the functions 
t^{G,v) and det A </3 (G, c) are never proportional if the genus o/E is greater 
or equal to two. 

Via the method developed in [8j Subsection 4.3], this theorem can be 
interpreted as a wide-reaching generalization of the main result of [61 [7] 
which was obtained via the Pfafhan method. In our opinion, this is a good 
example of how simpler and more natural a proof can get, when the Kac- 
Ward method is used instead of the Pfafhan one. (See Remark 14.71 below 
for a more detailed comparison.) In the case of the flat torus, the theorem 
above implies a relation between the free energy of the critical Z-invariant 
Ising model on G and the normalized determinant of the critical discrete 
Laplacian on G. The former quantity was computed by Baxter [3], the 
latter by Kenyon [27J, and our equality allows to obtain any of these two 
results MS M corollary of the other one. 

The paper is organized as follows. In Section [21 we define our model: the 
Ising model on graphs isoradially embedded in flat surfaces (Definition [1]), 
with critical weights (Definition [2]) . In Section O we introduce the </?- twisted 
Kac-Ward matrices for graphs in flat surfaces (Definition [3]) , and we show 
how they can be used to compute the Ising partition function (Theorem 13. 3|) . 
Section H] deals with the case of isoradially embedded graphs with critical 
weights, and contains our main results. We start with the combinatorial 
interpretation for the Kac-Ward determinants with critical weights (Propo- 
sition [4TT]). Then, we prove the equality relating the Kac-Ward determinants 
of dual isoradially embedded graphs (Theorem 14. 4|) . Also, we relate the Kac- 
Ward determinants with the determinant of the critical discrete Laplacian 
on G (Theorem I4.6[) . In a last paragraph, we explain how the Kac-Ward 
determinant can be understood as a discrete version of the <9-torsion of the 
underlying Riemann surface (Subsection 14. 4[) . 

Acknowledgments. This research was supported by the European Re- 
search Council AG CONFRA and by the Swiss NSF. Part of this paper 
was done at the ETH in Zurich, and it is a pleasure to thank the Math- 
ematics Department of the ETHZ for providing such an excellent working 
environment. The author also wishes to thank Cedric Boutillier and Hugo 
Duminil-Copin for comments on an earlier version of the manuscript, and 
Martin Loebl for valuable discussions. 

2. The critical Ising model on isoradial graphs 

The aim of this section is to explain the setup of the model that we will 
be studying in this paper: the critical Ising model on graphs isoradially em- 
bedded in a flat surface (Definitions [T] and [2]). To motivate this definition, 
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we start by recalling what is meant by high and low-temperature expan- 
sions for the Ising model, leading to the Kramers- Wannier duality argument 
(Subsection 12 . X [) . The models on which such an argument can be applied 
are called Z-invariant Ising models (Subsection 12. 2p . They are all defined on 
planar (or toric) graphs, but a generalization of these models to surfaces of 
arbitrary genus then naturally leads to our definition (Subsection 12. 3p . 

2.1. Kramers- Wannier duality. Let G be a finite graph with vertex set 
V(G) and edge set E(G). A spin configuration on G is a map a: V{G) — > 
{— 1, +1}. Any positive edge weight system J = (J e )eeE(G) on G determines 
a probability measure on the set £l(G) of such spin configurations by 

P ^ = Z jl (G) 6XP ( ^ JefT «^) ' 

where 

Z J {G) = ^2 exp ( ^2 JeQ-uO-^j 

o-en(G) e=(u,v)eE(G) 

is the partition function of the Ising model on G with coupling constants J. 
As observed by van der Waerden [40J, the identity 

exp( J e a u a v ) = cosh(J e )(l + tanh( J e )a u a v ) 
allows to express this partition function as 

Z\G) = ( J] cosh( J e )) J2 II i 1 + tanh( J e )a u cT,) 

eeE(G) aeU(G) e=(u,v)eE(G) 

= ( H cosh(J e ))2l y ( G )l J2 II tanh ( Je )' 

e€E(G) ' 7e£(G) eG7 

where £(G) denotes the set of even subgraphs of G, that is, the set of 
subgraphs 7 of G such that every vertex of G is adjacent to an even number 
of edges of 7. This is called the high-temperature expansion of the partition 
function. 

Let us now assume that the graph G is planar, and let G* denote its dual 
graph. For any spin configuration a 6 Q(G), consider the subgraph of G* 
given by all edges e* G E(G*) dual to e = (u, v) with a u / a v . Clearly, this 
is an even subgraph of G*. Furthermore, since G is planar, this defines a 
surjective map f2(G) — > £(G*) such that a and a' have same image if and 
only if a' = —a. This leads to the following low-temperature expansion of 
the Ising partition function: 

z\g) = 2 j2 n ex p( je ) n ex p(- 2j e ) 

7*e£(G*) eeE(G) e*£-y* 

=2( n ex P (j e )) y, n ex p(- 2j e )- 

eGE(G) 7*g£(G*)e*e7* 
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Thus, if we assign weights J to E{G) and J* to E{G*) in such a way that 
tanh(J e *) = exp(— 2J e ), or more symmetrically, 

sinh(2J e )sinh(2J e *) = 1, 

we obtain that the partition functions Z J (G) and Z J * (G*) are proportional 
to each other. 

In the case of the square lattice with constant weight system J, this is 
enough to determine the critical value J c of the coupling constant. Indeed, 
assuming that the free energy of the model is analytic everywhere except at 
a single point, this point must be equal to J c and to J*, since the square 
lattice is self-dual. The equality J c = J* then leads to the explicit value 

J c = log Vl + V2. 

2.2. The Z-invariant Ising models. This beautifully simple argument, 
the celebrated Kramers- Wannier duality [31J , is not sufficient to determine 
the critical value of the coupling constant on a graph that is not self-dual: 
it only relates this critical value to the one for the dual graph. However, for 
some planar graphs, this can be obtained "with little additional labor" |41j . 

Let us start with the example of the hexagonal lattice H, and let H' 
denote the graph obtained from H by a star-triangle transformation at a 
vertex v as illustrated in Figure [TJ Assume one can assign coupling con- 
stants J to the edges of H and J' to the newly created edges of H' in such 
a way that Z J '(H') = R v Z J {H) for some function R v of the coupling con- 
stants of the edges around v. Since the triangular lattice T can be obtained 
from H by such transformations, it would follow that Z K (T) = RZ J (H) for 
some controlled R and well-chosen coupling constants K, J. This, together 
with the Kramers- Wannier duality, would lead to an equality of the form 
Z K (T) = kZ K *(T), with K 1— > K* some involution and k an explicit func- 
tion of K and K*. Arguing as above, the critical points should be self-dual 
under this involution, leading to the equality k = 1 and an exact description 
of these critical points for the triangular and hexagonal lattices. 

This strategy of using invariance under the star-triangle transformation 
(or, Z-invariance) can be applied not only to the hexagonal-triangular lat- 
tices, but to a wide class of planar graphs [3 0]. It turns out that this class 
of graphs on which a Z-invariant Ising model can be defined coincides with 
the graphs that admit an isoradial embedding in the plane lllj: this is 
an embedding such that each face is inscribed in a circle of radius one, with 
the circumcenter in the closure of the face. Furthermore, the corresponding 
critical coupling constants admit a very simple geometric description: they 
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Figure 2. An edge e of an isoradial graph, the associated 
rhombus, and the half-rhombus angle 9 e . 



are given by 



1 , / 1 + sin I 
J e = - log 



2 \ cos 9 e 

where 9 e 6 (0, vr/2) is the half-rhombus angle associated to the edge e, 
as illustrated in Figure [2J For example, the square lattice is isoradially 
embedded with all half-rhombus angles equal to 9 = 7r/4, leading to the 

critical coupling constant J c = log yl + \/2 as above. On the other hand, 
the triangular and hexagonal lattices are isoradially embedded with angles 
9 = tt/6 (resp. vr/3), so the corresponding critical values are equal to J c = 



logvv3 (resp. log y2 + \/3). 

This geometric description of the critical coupling constants becomes even 
nicer when using the high-temperature expansion, as we obtain the critical 
weights 

Wt n , /l. / 1 + sin 9 e \ \ 1 — cos6* e . .. 

u e = tanh( J e ) = tanh - log — = — = tan(0 e /2). 

V 2 V cos dp / / sm 6L 



2.3. Isoradial graphs in flat surfaces. In the present paper, we will study 
more general models where the graph G is not assumed to be planar. The 
proper generalization of planar isoradiality is obtained by considering so- 
called flat surfaces with cone-type singularities. Let us quickly recall their 
definition and main properties, referring to |39j for further details. 

Given a positive real number ■d, the space 

C$ = {(r,t) : r > 0, t G R/0Z}/(O,t) ~ (0,0 

endowed with the metric ds 2 = dr 2 + r 2 dt 2 is called the standard cone 
of angle Note that the cone without its tip is locally isometric to the 
Euclidean plane. Let S be a surface with a discrete subset S. A flat metric 
on S with cone-type singularities of angles {"& x }x£S supported at S is an atlas 
{<f>x '■ U x — > U' x C Cfi^xss, where U x is an open neighborhood of x £ S, cj) x 
maps x to the tip of the cone C$ x , and the transition maps are Euclidean 
isometries. 

This seemingly technical definition should not hide the fact that these ob- 
jects are extremely simple and natural: any such flat surface can be obtained 
by gluing polygons embedded in R 2 along pairs of sides of equal length. For 
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example, a rectangle with opposite sides identified will define a flat torus 
with no singularity. On the other hand, a regular 4g-gon with opposite sides 
identified gives a flat surface of genus g with a single singularity of angle 
2ir{2g — 1). In general, the topology of the surface is related to the cone 
angles by the following Gauss-Bonnet Formula: if E is a closed flat surface 
with cone angles {$ x }x<=s, then 



where x(X) is the Euler characteristic of X. 

Definition 1. A graph G is isoradially embedded in a flat surface X if the 
following conditions are satisfied: 

- X is a compact orientable flat surface with cone-type singularities; 

- each edge of G is a straight line in X; 

- each closed face / of G C X contains an element Xf at distance 1 
from all vertices of df; 

- a singularity of X is either a vertex of G or a vertex x f of the dual 
graph G*, that is, the singular set is contained in V{G) U V(G*). 

Given an isoradially embedded graph G C X, each edge e 6 E(G) has 
an associated rhombus as illustrated in Figure [2j Therefore, the metric 
space X should simply be understood as rhombi pasted together along their 
boundary edges. This observation also leads to the following fact. 

Proposition 2.1. Any finite graph G can be isoradially embedded in aflat 
surface. 

Proof. Fix an arbitrary angle 8 S (0,7r/2) and associate to each half-edge e 
of G the isosceles triangle illustrated below. 



For each vertex v G V(G), choose a cyclic ordering of the d v half-edges 
of G coming out of v and build the associated star St(y) as follows: glue 
together the corresponding d v triangles following the chosen cyclic ordering 
around v. Note that St(v) is a flat surface with one singularity of cone angle 
29d v , and that the cyclic ordering endows this surface with an orientation. 
For each edge e = (v,w) £ E(G), glue together the stars St(v) and St(w) 
along their boundary by pasting together the two triangles associated to e 
in the unique way that is consistent with the orientation of the stars. (If e 
is a loop at v, just glue together the corresponding sides of the star St(v).) 
By construction, G is isoradially embedded in the resulting metric space X, 
which is a compact oriented flat surface. □ 



£(2ir - 0*) = 27r X (E), 



ices 




Following the discussion of the previous subsections, we shall adopt the 
following terminology. 
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Definition 2. Let G be a graph isoradially embedded in a flat surface E. 
The critical weight associated to the edge e G E(G) is defined by 



where 9 e G (0, 7r/2) is the half-rhombus angle associated to the edge e. The 
partition function for the critical Ising model on G is given by 



where £(G) denotes the set of even subgraphs of G. 

3. The Kac-Ward formula for graphs in flat surfaces 

In jS], we gave a generalized Kac-Ward formula for the Ising partition 
function on any finite weighted graph (G, x). The aim of this section is to 
show that when the graph G is embedded in a flat surface, the generalized 
Kac-Ward matrices take a particularly simple form - whatever the weight 
system x on G is. In the next section, we will consider the case of isoradial 
graphs with critical weights. 

3.1. Kac-Ward matrices for graphs in flat surfaces. Let us start with 
some general terminology and notation. Given a weighted graph (G, x), let 
E = E(G) be the set of oriented edges of G. Following [37], we shall denote 
by o(e) the origin of an oriented edge e E E, by t{e) its terminus, and by e 
the same edge with the opposite orientation. By abuse of notation, we shall 
write x e = x g for the weight associated to the unoriented edge corresponding 
to e and e. 

Now, assume that G is embedded in an orientable flat surface S so that 
each edge of G is a straight line, £ \ G consists of topological discs, and the 
set S of cone-type singularities is contained in V{G) U V(G*). As above, let 
$ x denote the cone angle of the singularity x G S. Fix a unitary character 
ip of the fundamental group of S, that is, an element of 



Hom(vri(S), C/(l)) = Hom(7r 1 (S), S 1 ) = Hom(#i(E), S 1 ) = H^^S 1 ). 



As G C S induces a cellular decomposition of E, one can represent such 
a cohomology class by a cellular 1-cocycle, that we shall also denote by (p. 
This is nothing but a map from the set E of oriented edges of G into S , such 



where /3(e',e) G (0, denotes the angle from e' to e, as illustrated in the 
left part of FigureO We shall call the matrix I—T v the tp-twisted Kac- Ward 
matrix associated to the weighted graph (G, x), and denote its determinant 



v e = tan(# e /2) 



s(g»= e n 



7G£(G) e€7 




by tV(G,x) = det(J-T^). 
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Figure 3. The angles /3(e',e) and a(e,e') = it - /3(e',e). 



Remark 3.1. If t{e) = o(e') = v is not a singularity (that is, if the cone angle 
•d v is equal to 2ir), then the complex number iexp (— |/3(e',e)) is equal to 
exp (|a(e, e')), with a(e,e') G (— vr,7r) the angle from e to e' illustrated 
in the right part of Figure O If v does belong to S, then the dotted line 
drawn there does not make sense anymore, hence the necessity to adopt this 
slightly less intuitive definition. 

3.2. Discrete spin structures on flat surfaces. The Kac-Ward matrices 
are particularly useful when ip is a special type of 1-cocycle, namely a discrete 
spin structure. In this paragraph, we shall recall the definition and main 
properties of these objects, slightly generalizing Section 3.1 of [9]. 

Loosely speaking, a spin structure on an oriented surface E is a way 
to count parity of rotation numbers for closed curves in E. In the plane, 
there is a unique way to do so, and therefore a unique spin structure. On 
a closed orientable surface of genus g however, there are exactly 2 2g dis- 
tinct spin structures. More precisely, one can identify the set S(E) of spin 
structures on E with the set of quadratic forms on E [23]: these are the 
maps q: i?i(E;Z2) — > Z2 such that q(x + y) = q(x) + q(y) + x ■ y for all 
x,y G i?i(E;Z2), where x ■ y denotes the intersection number of x and y. 
Note that the difference of two quadratic forms is a linear form. Therefore, 
this set admits a freely transitive action of the abelian group i/ 1 (E; Z2); in 
other words, it is an affine -ff 1 (E; Z2)-space. The Arf invariant of a spin 
structure is defined as the Arf invariant of the associated quadratic form, 
that is, the number Arf (5) G Z2 satisfying 

(_l)Arf(g) _ J_ ^ (-1)90). 

Coming back to flat surfaces, let us assume that G is a graph embedded 
in a flat surface E so that E \ G consists of topological discs, and let X de- 
note the induced cellular decomposition. We shall now explain how, in such 
a situation, it is possible to encode spin structures on E by some cocycles 
A G Z 1 (X; S 1 ). Let us assume that all cone angles $ x are positive multiples 
of 2-7T, i.e. that E has trivial local holonomy. Then, the holonomy defines 
an element Hoi of Hom^E), S 1 ) = il^E; S 1 ) = H^X; S 1 ). We shall call 
a cocycle k G Z 1 (X;S 1 ) such that [k™ 1 ] = Hoi a discrete canonical bundle 
over E. Note that such a cocycle k is very easy to determine. Indeed, it is 
always possible to represent E as planar polygons P with boundary identi- 
fications. Furthermore, these polygons can be chosen so that G intersects 
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dP transversally, except at possible singularities in Sfl V(G). Define k by 



where 9 denotes the angle between the sides of dP C C met by the edge 
e. If S n V{G) is empty, this defines completely a natural choice of dis- 
crete canonical bundle k. Otherwise, the partially defined k above can be 
extended to a cocycle yielding a discrete canonical bundle. 

Mimicking the continuous case (in the version developed by Atiyah [2]), 
let us define a discrete spin structure on S as any cellular 1-cocycle A £ 
Z 1 (X;S 1 ) such that A 2 = k. Two discrete spin structures will be called 
equivalent if they are cohomologous. The set §(X) of equivalent classes of 
discrete spin structures on S is then given by 



Note that if the flat surface £ has trivial holonomy, then [«] is trivial, so the 
set §(X) is equal to the 2g-dimensional vector space i? 1 (£;Z2). In general, 
GcE can be described via planar polygons as explained earlier. In such a 
case, and assuming that the singular set S n V(G) is empty, a discrete spin 
structure is given by 



where exp(— iO/2) denotes one of the square roots of the angle between the 
sides of dP C C met by the edge e. 

One easily checks that the set S(X) is an affine Z2)-space. Fur- 

thermore: 

Proposition 3.2. If all cone angles o/S are odd multiples of2ir, then there 
exists a canonical H 1 (T>;Z2)-equivariant bijection §(X) — > §(£)• 

Proof. Let k £ Z 1 (X;S 1 ) be a fixed discrete canonical bundle over S. For 
each A £ Z (X\ S 1 ) such that A 2 = k, we shall now construct a vector field 
V\ on £ with zeroes of even index. Such a vector field is well-known to define 
a spin structure, or equivalently - by Johnson's theorem [23] - a quadratic 
form q\ on Z2). The proof will be completed with the verification that 

two equivalent A's induce identical quadratic forms, and that the assignment 
[A] 1 y q\ is Z2)-equi variant. 



First replace the cellular decomposition X of E by X', where each singu- 
larity v £ V(G) n S is removed as illustrated above. Obviously, A induces 
A' £ Z 1 (X';S' 1 ) by setting A'(e) = 1 for each newly created edge e. Fix 
an arbitrary orientation u) of the edges of X' . This allows to represent 
A' £ Z l (X';S 1 ) as follows: write A'(e) = exp(i/3 A /(e)) with < /3y(e) < 2vr 




1 if e is contained in the interior of P; 

exp(—i9) if e meets dP transversally, 



§(X) = {[X]eH 1 (X;S 1 )\[\] 2 = [ K ]}. 




1 if e is contained in the interior of P; 

exp(—i9/2) if e meets dP, 
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if e is an edge oriented by u), and set /3y(e) = —(3\>(e) for the reverse edge. 
Fix an arbitrary tangent vector V\(v) at some arbitrary vertex v of X', and 
extend it to the 1-skeleton G' of X' as follows: running along an edge e, 
rotate the tangent vector by an angle of 2/3y(e) in the negative direction. 
Since A is a cocycle such that A 2 = k, and since each cone angle is a mul- 
tiple of 2tt, this gives a well-defined vector field along G' . Extend it to the 
whole surface X by the cone construction, creating one zero in each face of 
X and at each element of V{G) n S. Obviously, the resulting vector field 
V\ depends on the choice of uj, but not in a crucial way. Indeed, reversing 
the orientation ui on a given edge e either does nothing (if A'(e) = 1), or 
corresponds to adding two full twists to the vector field along e. Therefore, 
the parity of winding numbers with respect V\ are independent of oj. In 
particular, one easily checks that a zero of V\ is of even index if and only if 
the corresponding cone angle is an odd multiple of 2ir, which we assumed. 

As explained in [24J, the quadratic form q\ : fli(£; Z2) — > 1>2 correspond- 
ing to V\ is determined as follows: for any regular oriented simple closed 
curve C C £ \ S, the number </a([C]) + 1 is equal to the mod 2 winding 
number of the tangential vector field along C with respect to the vector 
field V\. For an oriented simple closed curve C C G, we obtain the following 
equality modulo 2: 

q x ([C}) = 1 + ^(E &(e) + «^)) • 

eCC veC 

where the first sum is over all oriented edges in the oriented curve C, and 
a v (C) is the angle illustrated below. (This angle should be interpreted 
as explained in Remark 13.11 ) Obviously, equivalent A's induce the same 




quadratic form q\. Finally, given two discrete spin structures Ai,A2, the 
cohomology class of the 1-cocycle A1/A2 £ Z 1 (X; {±1}) is determined by its 
value on oriented simple closed curves in G. For such a curve C, we have 

(Ai/A 2 )(C) = exp (^(^(e) - /3 A2 (e))) = exp (ivr(g Al - </a 2 )([C])). 

eCC 

Therefore, the assignment [A] 1— > q\ is H 1 ^; Z^-equivariant, which con- 
cludes the proof. □ 



3.3. The Kac-Ward formula for flat surfaces. We are finally ready to 
prove the main result of this section, motivating the introduction of twisted 
Kac-Ward matrices and discrete spin structures. As before, let (G,x) be a 
weighted graph embedded in a closed orientable flat surface E so that each 
edge of G is a straight line, S \ G consists of topological discs, and the set 
S of cone-type singularities is contained in V(G) U V(G*). 
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Theorem 3.3. If all cone angles are odd multiple of2ir, then the Kac-Ward 
determinant r ¥ '(G, x) is the square of a polynomial in the x e 's whenever p> 
is a discrete spin structure on G CS. Furthermore, if t v {G,x) 1 ^ 2 denotes 
the square root with constant coefficient equal to +1, then the Ising partition 
function on G is given by 

z ^*) = y 9 E (-i^Mg,*) 1 /*, 

Aes(x) 

where g is the genus of S and Arf (A) 6 ^2 the Arf invariant of the spin 
structure corresponding to A. 

The demonstration given below is by no means self-contained: it only 
consists in recasting the flat surface case in the more general (and more 
complicated) topological setting discussed by the author in [8j. We refer to 
this article for further details. 

Proof. By Bass' Theorem [3], r^G, a?) is given by 

T *(G,x)=det(/-T»')= n i 1 - II T ^')> 

7 eIP(G) (e,e')£7 

where 7{G) is the (infinite) set of prime reduced oriented closed paths in G, 
and the second product is over all pairs of consecutive oriented edges in the 
oriented path 7. By definition of T^, 

II T e,e> = <p(rr) ex P (l a (7)) a?(7)j 

(e,e')67 

where ^(7) = n e g7 X e and 0(7) is the sum of the angles a(e,e') along 7 
(interpreted as in Remark 13 . X [> . This equality already shows that t v (G,x) 
does not depend on the choice of the 1-cocycle representing the cohomology 
class (p. Furthermore, if </? G Z 1 (X;S 1 ) is a discrete spin structure, then 

(93(7) exp (|a(7))) 2 = vdf exp(ia(7)) = k' 1 ^) = 1 

for all oriented closed path 7 in G. This implies that ^(7) exp (f 01(7)) 
always belongs to {±1}- Therefore, for such a ip = X £ §>(X), 

t\g, X )= n (i-(-ir h) x(7)) 2 , 

for some u^a(t) £ where the equivalence relation on 5 (G) is given by 
7 ~ —7. Using the notations of the proof of Proposition 13.21 the element 
wx("f) satisfies 

(_1)«a(7) = A( 7 ) exp (§a( 7 )) = exp (*X>(e) + \ EM7))- 

eC7 U67 

Therefore, Wa(7) is nothing but the mod 2 winding number of the tangent 
vector field along 7 with respect to the vector field V\ associated to the 
discrete spin structure A. The formula now follows from Proposition 13.21 and 
[SI Corollary 2.2]. □ 
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Remark 3.4. The whole setting can be extended to encompass graphs em- 
bedded in flat surfaces with boundary. If there is exactly one boundary 
component, then Theorem 13.31 extends verbatim. (In particular, it applies 
to domains in the plane, where this theorem is exactly the original Kac- 
Ward formula [31].) If the flat surface has several boundary components, 
then the formula is slightly more complicated. 

4. The Kac-Ward matrices with critical weights 

As proved in the previous section, the ^-twisted Kac-Ward matrices can 
be used to compute the Ising partition function for any weighted graph em- 
bedded in a flat surface. We shall now assume the graph to be isoradial and 
the weights to be critical (recall Definitions [JJ and . We will start with 
a combinatorial interpretation for t^(G,v) (Proposition 14. that we then 
use for two of our main results. First, we prove a duality theorem relating 
tV{G,v) and t v (G*,v*) (Theorem H3D . Then, we show that t^(G,v) co- 
incides up to a multiplicative constant with the determinant of the critical 
discrete Laplacian on G if and only if the genus of E is zero or one (Theo- 
rem |T6|). In a last subsection, we explain how t v {G,v) can be understood 
as a discrete version of the 9-torsion of the underlying Riemann surface. 

4.1. A combinatorial interpretation for t v (G, u). Let us start with 
some notations. 

Given a graph G, let 3~(G) denote the set of subgraphs F C G such that 
F spans all vertices of G, and no connected component of F is a tree. Also, 
for a graph F embedded in an oriented surface E, we shall denote by N(F) 
a small tubular neighborhood of F in S. (This is simply F "thickened", as 
illustrated in the middle of Figured!) Since E is oriented, so is N(F), and 
this induces an orientation on the boundary dN(F) of N(F). Therefore, 
dN(F) consists of a disjoint union of oriented simple closed curves 7 on E. 
This is illustrated in Figured) 

Proposition 4.1. Let G be a graph isoradially embedded in a flat surface 
E, and let us assume that all cone angles i9 v of singularities v £ V(G) are 
odd multiple of2n. Let v e = tan(# e /2) denote the critical weight system on 
G C E, and set fi e = itan(0 e ). Then for any 1-cocycle (p, 

t*(G,v) = C II (l-p(7))/*(n 

F&{G) -yCdN(F) 

where the product is over all connected components 7 of the ( clockwise ori- 
ented) boundary of a tubular neighborhood N(F) of F in E, n(F) = Y\ eeF fJ, e , 
and the constant C is equal to 

c=<-u»w> n «**./«> n jf^m- 

with x(G) = \V(G)\ - \E(G)\ the Euler characteristic of G. 

The proof being quite substantial, we will split it into two lemmas. Let 
us begin with some notation. We shall write V = V(G) and E = E{G) 
for the sets of vertices and edges of the graph G, E = E(G) for the set 
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Figure 4. Part of a graph F embedded in a surface S, the 
corresponding tubular neighborhood N(F), and the oriented 
boundary dN(F). 



of oriented edges of G and £(E) for the complex vector space spanned by 
E. Obviously, E can be partitioned into E = \_\ v& yE v , where E v contains 
all oriented edges e with origin o(e) = v. Now, let us cyclically order the 
elements of E v by turning counterclockwise around v. (As £ is orientable, 
this can be done in a consistent way.) Given e £ E v , let R(e) denote the 
next edge with respect to this cyclic order, as illustrated in Figure This 
induces an endomorphism R of £(E). Also, let J denote the endomorphism 
of £(E) given by J(e) = e. Finally, we shall write /i for the endomorphism 
of £(E) given by /i(e) = [i e e, and similarly for any weight system and for <p. 

The aim of the first lemma is to relate r ¥ '(G, u) to the determinant of a 
more tractable matrix. 



Lemma 4.2. For G C S as in Proposition \4-l\ and for any 1-cocycle cp, 

with M = I + Jipfi - R{fj, + 1) G End(L(E)) and X (G) = \ V(G)\ -\E(G)\. 

Proof. Let Succ £ End(L(K)) be defined as follows: if e is an oriented edge 
with terminus t(e) = v, then 

Succ(e) = iip{e)v e u>(e',e)e', 

e'eE v 

where u)(e', e) = exp(— |/3(e', e)) for e' ^ e £ E v with /3(e', e) as in Figure[3l 
and w(e, e) = —1. Also, let S £ End(L(E,)) be the endomorphism given 
by S = Succ + iJ(pv. By definition, the y-twisted Kac-Ward matrix is the 
transposed of / — S. Now, consider the matrix 

A = (I — S)(I + iJipu) = I — Succ + Com, 

where 

Com(e) = —i(p(e)v e S(e) = u(e',e)e' 

e'eB v 

if e has origin o(e) = t> . Since 

de t( , + _ n det ( 1 ^ ) . na + ,1), 
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Figure 5. The endomorphism R (to the left), and the equal- 
ity f3(e, R(e)) =9 e + 6 R{e) (to the right). 



we get the equality 

(1) T*(G,v) = H(l + v 2 e r 1 detA. 

(The computation above is a variation on a trick due to Foata and Zeilberger, 
see [HI Section 8].) 

Let Q G End(L(E)) be defined by 

Q(e)=exp{ip(e,R(e))) e, 

and set P = I — RQ. Obviously, this endomorphism decomposes into P = 
(& veV P v with P v £ End(L(E v )), and one easily computes 

det P v = 1 - Y[ ex P (20(e, R(e))) = 1 - exp (%& v ) = 2, 

e&E v 

as the cone angle $ v is an odd multiple of 2ir. Hence, the determinant of 
P is 2^1 The point of introducing this P is that it can be used to greatly 
simplify the matrix A. Indeed, let us compute the composition PA. If e has 
terminus v, then 

PSucc(e) = (I — RQ)i(p(e)u e u(e',e)e' 

e'eE v 

= i<p(e)v e (w(e , ,e)-o;( J R- 1 (e / ),e)exp(§/3(i?- 1 (e , ),e')))e' 

e'€E v 

= -2iip(e)v e (e). 

(To check that the coefficient of e' = R(e) vanishes, we use once again the 
fact that is an odd multiple of 2n.) Therefore, we have the equality 
PSucc = —2iJ(pv. Similarly, given e with o(e) = v, 

PCom(e) =(J - RQ)v 2 e ^ u}(e',e)e' 

e'€E v \{e} 

=vl M e '' e ) " e) exp (§/3(irV), e'))) e' 

e'eB v \{e,ie(e)} 

+ ^ e 2 (w(ii(e),e) J2(e) - w(i2 _1 (e),e) exp (i/3(i2 _1 (e), e)) e) 
= -^ 2 (/ + J RQ)(e). 
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These two equalities lead to 

PA = P(I - Succ + Com) 

= (J — RQ) + 2iJipv -{1 + RQ)v 2 
= (1 - v 2 )I + 2%Jipv - RQ(1 + v 2 ). 

We finally get 

(2) det A = 2"l' l/ l Y[(l-v 2 ) 2 detA', 

where A' is given by 

A' = I + iJip -^-s - RQ 

1 — u z 1 — v z 

It is now time to use the fact that the weights v e are not any weights, but 
the critical ones given by v e = tan(0 e /2). First observe that 



= itan{O e ) = fx e and 



1 — cos{a e ) 

Next, note the equality /3(e, R{e)) = 9 e + #_R( e ) illustrated to the right of 
Figure [5j This implies that RQ(e) = exp (§0 e ) exp (§0R( e )) -R(e)- Multi- 
plying each column of A' (corresponding to e) by exp (§# e ) and each line 
(corresponding to e') by exp (— |0 e '), we obtain a new matrix M with 

(3) detA' = detM and M = J + J<^ - i?(l + fx), 

since # g = # e and e ^^ e ^ = 1 + fx e . Equations ([TJ, ^ and © give the 
statement of the lemma. □ 

In a second lemma, we now give a combinatorial interpretation of the 
determinant of this matrix M. 

Lemma 4.3. The endomorphism M = 1+ Jipfx — R(l + fx) o/£(E) satisfies 
detM = (-l)M + Me) £ II (1 " P(7)) M(n 

eG-E KG 1 cdN(F) 

V(F)=V x ' 

the sum being on all subgraphs F of G spanning all vertices of G. 

Proof. By definition, the coefficients of M are given by 

\ ife = e'; 

M ej£ / = < (p(e)fx e if e' = e; 

k -(l + A*e) Xe? = R(e), 
and vanish otherwise. Let us compute directly the determinant of M as 

detM= £ (-l) Sgn(CT) II M ^)- 

o-eS(E) eGE 
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Each permutation a G 5(E) decomposes into disjoint cycles, inducing a 
partition E = | |^ Ej(a) into orbits of length £j(a). The corresponding con- 
tribution to the determinant is 

( _i).«nw n^(e)=n(- i )^ (<7)+i n 

Since M e ^ e is equal to 1, the oriented edges that are fixed by a contribute a 
trivial factor 1 to this product and the corresponding orbits can be removed. 
By definition of M, a permutation a 6 5(E) will have a non-zero contribu- 
tion only if each of the n(a) remaining orbits forms a cycle of oriented edges 
of G such that each oriented edge e is either followed by J(e) = e or by 
R(e), and such that these cycles pass through each edge of G at most twice, 
and if so, in opposite directions. If T(G) denotes the set of such union of 
cycles, we have 

detM= y, (-ir (7) nn- M ^(e). 

7er(G) j ee 7j 

where each 7 G r(G) is written as a union of cycles 7 = (Jji7i lj- Now, 
for any fixed 7, a given edge e e E will fall in one of the following five 
categories: 

(i) e is covered by 7 in both directions, as part of a cycle of the form 
(. . . , R~ 1 (e), e, e, R(e), . . . ); the corresponding contribution to the 
determinant of M is — (p(e)fj, e (l + fx e ). 

(ii) e is covered in both directions by the cycle (e, e), so the contribution 
is — /Ug. (The minus sign comes from the contribution of this cycle 
to 71(7).) 

(m) e is covered by 7 in both directions, as part of cycles of the form 
(. . . , i? _1 (e), e, R(e), . . . , i? _1 (e), e, R(e), . . . ); in this case, the con- 
tribution is (1 + /x e ) 2 . 
(ivj e is only covered in one direction, so the contribution is (1 + fi e ). 
(v) e is not covered at all, and the contribution is 1. 

For any given element 7 £ T(G), any edge of type (ii) in 7 can be removed 
(i.e. replaced by an edge of type (v)) and the resulting union of cycles 
will still belong to T(G). The converse also holds: any type (v) edge in an 
element of T(G) can be replaced by a type (ii) edge, the result will be in 
T(G). Therefore, each time an edge appears as an edge of type (ii) of some 
element of T(G), it also appears as an edge of type (v) of some other element 
of T(G) and vice versa. Using the equality (1 — y? e ) = (1 + /x e )(l — ;U e ), we 
therefore can factor out a term 1 + /j, e for each e££, leading to 

detM=na+^) £ (-ir {7) n fW-ve) n ^ n 

e£E 7er(G) e€7(t) eey(ii) eey(iii) 

where j(i) (resp. j(ii) , j(iii)) denotes the set edges of type (i) (resp. 
(ii), (Hi)) of 7. Next, we wish to expand the last product above as 

n (i+^e)= £ n^- 

eS7(m) 7'C7(m) eg7' 



The critical Ising model via Kac-Ward matrices 



19 



Using the notation 

f(G) = {( 7 ,V)l7er(G), 7 'C7(m)}, 
we get the equality 

detM=n(l+^e) E (-1)^V(7) II Me, 
e&E 7G?(G) eeSupp<a) 

where 71(7) = 72(7) + |7(i)|, the support of 7 is Suppi^y) = 7(7) U 7(ii) U 7' 
and 99(7) = n e G7(i) ^( e ) f° r an Y 7 £ r(G). In other words, 

(4) - d * M <M = E c FM(n where 4= £ (-1)^(7). 

We shall now check that Cp vanishes whenever F does not span all vertices 
of G. Let us assume that the vertex v is not spanned by F. Let 7^ be 
the cycle (e^i?(e), i? 2 (e), . . . , i2 d_1 (e)) with e £ E v and d = |£7„|. The 
set of 7 G T(G) with Supply) = F is equal to the disjoint union of To 
and Ti, where To (resp. Ti) denotes the set of such 7 = (7,7') with 7 
containing (resp. not containing) j v . As v does not belong to Supply), the 
mapping (7,7') 1— > (7 U 7^,7') gives a well-defined bijection /: To — > T\. 
Since 71(7(7)) = 71(7) + 1 and <p(f(j)) = <p(j), we get 

4 = E (-!) n( ^V(7) = E (c- 1 )^ + (-ir (/W) M7) = 0. 

7er uri 7er 

So, let us assume that F C E spans all vertices of G, and let {T/Jfc denote 
the connected components of F viewed as a subgraph of G. Given 7 G 
let 7I denote the restriction of the cycles composing 7 to the oriented edges 
of 2V One easily checks that if (7,7') G has support equal to F, then 
7|T fc belongs to T(T k ) and 7|T fc (») = 7(*) HTj for • G Hence, 
sending (7,7') to (7^, 7' n gives well-defined maps 

7T fc : {7 G F(G) I 5«pp(7) = F} -> {7 fc G f (T fc ) | Supp{%) = T k }, 

which in turn induce 

{7 G f(G) I Supp(j) =F}^ Y[tf k G F(r fc ) I Suppftk) = 

Using the fact that F spans all vertices of G, one can check that this map 
is a bijection. Since n(j) = ^2 k n(TTk(j)) and 92(7) = JIa: V'( 7r fc(7)) f° r an Y 7 
with support equal to F, it follows that 

(5) 4 = U c n where c£ = £ (-1)^^(7). 

fc 7er(T) 

5upp(7)=T 

Note that the value of Cj, does not depend on G anymore. Furthermore, 
since Supply) must be equal to the whole of T, 7 = (7,7') must satisfy 
7' = 7(m). Therefore, 

4= E (-i) B(7)+l7(i)l v(7), 

7er(T) 

Supp(i)=T 
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where Supp(^) = 7(4) U U 7 (Hi) and ^(7) = n e e 7 (i) ^( e )- 

Let iVy G £W(£(E(T)) be the endomorphism given by Nj, = R — Jip. 
Developing the determinant of Nj, explicitely, we get 

detJV£ = Yl (-l) n(7)+l7(i)l ^(7) = 4. 
7er(T) 

Supp(~f)=T 

As before, the set E(T) is equal to |_|i>gV(T) E V (T), and the endomorphism 
R splits into R = 0„ e ^ (T ) R v - Since deti^ = (-1)I^( T )I+1, we get 

4 = det(-R) det(/ - R~ l J<p) = (-1)1^)1 det(J - R^Jip). 

Now, the orbits Ej of the action of R~ 1 J on the set E(T) correspond ex- 
actly to the connected components of the oriented boundary of a tubular 
neighborhood N(T) of T in S, where N(T) is endowed with the clockwise 
orientation. Hence 

(6) 4=(-i)"^"n(i-n*'( e ))=(- i ) |vC01 n a-^)- 

j eeBj ' 7 CdiV(T) 

Since |F| = = Ek\ V ( T k)\ and JV(f) = Uk N ( T k)> Equations ©, 

© and (|6j) give the statement of the lemma. □ 

Proof of Proposition \4-l\ Lemmas 14.21 and 14.31 give the equality 

t*{G,v) = C II (l-¥>(7)Mn 



where the constant C is equal to 



COS( 



This is equal to the value given in the statement of the proposition, since 
II exp(^ e ) = II II ex P(^/2) = II exp(*V4). 

e£-B v€V e€E v v£V 

Finally, note that if a component of F is a tree, then dN(F) will contain 
a trivial cycle and the corresponding coefficient will vanish. Therefore, we 
can sum over all elements of the set 3~(G) defined at the beginning of the 
section. This concludes the proof of the proposition. □ 

4.2. The duality theorem. We are now ready to state and prove one of 
our main results. 

Theorem 4.4. Let G be a graph isoradially embedded in a flat surface S, 
and let v be the critical weight system on G. The dual graph G* is also 
isoradially embedded in E, and therefore admits a critical weight system v* . 
If all cone angles are odd multiples of '2tt, then for any ip G H 1 ^; S 1 ), 

2 \V(G*)\ Yl (l+cos(9 e *))r^(G*,u*) = 2l y ( G )l Yl (l+cos(0 e )) r^(G». 

e*eE(G*) eeE(G) 



The critical Ising model via Kac-Ward matrices 

1 



21 



Figure 6. A graph F (in solid lines) together with the asso- 
ciated graph vp(F) (in dashed lines), and the corresponding 
decomposition £ = N(F) U N(ip(F)). 



Proof. Given a subgraph F of G, let ip(F) denote the subgraph of G* given 
by 

xj)(F) = {e* e E(G*)\e(^ F}. 

Obviously, ifi defines a bijection from the set of subgraphs of G onto the set 
of subgraphs of G*, with inverse ^(F*) = {e£ E(G) \ e* <£ F*}. Let 7(G) 
denote the set of subgraphs of G spanning all vertices of G, and containing 
no cycle that is the boundary of a face. These two conditions being dual to 
each other, the map tp defines a bijection tp: 3~(G) — > 3"(G*). Now, consider 
the sum 



with c F 



F&{G) 



n (i 

~/CdN(F) 



Since dual edges have rhombus half-angles related by 

Therefore, 

= E ( 

F*£?(G*) 



S - 6» the 



weight fj, e = itan(# e ) satisfies \x e * = —fi e 



■1)1^)14 //(F*). 



eg£(G) 



Fe?(G) 



Furthermore, for any F £ 3~(G), one can decompose the surface £ as the 
union of two tubular neighborhoods N{F) and N(ip(F)) of i 7 and ip(F), 
pasted along their common boundary. (This is illutrated in Figure El ) In 
particular, the oriented boundaries satisfy dN(F) = —dN(il>(F)). There- 
fore, since ip is a 1-cocycle, 



: II (i-rt-y*)) 

7*ca/V(F*) 

: (_-q|£W(F)| -q 



n i 1 

7C9AT(F) 



^(7)4 = (-i)i° Jv ^i v - 1 (flj\r(F)) 4 



7C9Af(F) 



1) 



\8N(F)\ <p 



where \dN(F)\ denotes the number of boundary components of N(F). Note 
also that, since F spans all vertices of G, adding or removing an edge to F 
changes the parity of \E(F*)\ and of \dN(F)\. Therefore, the parity of their 
sum does not depend on F, and it is easily seen to be equal to the parity of 
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|V(Gr*)|, the number of faces of G C X. We have proved: 

]^[ /j,J 1 D' p (G, //) = (-l)^ G ^D^(G*,n*). 

e£E(G) 

Since Cp vanishes whenever F G 3~(G) does not belong to 3~(G), Proposi- 
tion 14.11 gives the equalities 

T =(-i)i^)i 2 -*( G ) jj exp( ^ /4) jj rTcSn^^'^ 

ueV(G) eeE(G) COS ^ e ' 

r * =( _ 1) |V( G -)| 2 -x(G-) JJ exp( ^ /4 ) J] cos(ge«) jyp^^ 

±A - ±A - 1 + COS t'»» 

dGV(G*) e*e£(G*) v ' 

where r and r* stand for r^G, i/) and r l|9 (G*,^*), respectively. The three 
equations above yield the equality 

2 \V(G*)\ JJ (l + cos(0 e .))r* = ^.2l y ( G )l JJ (l + cos(0 e ))r, 

e*sB(G*) ee£(G) 

where the constant A is equal to 

V ; n^exp(z^*/4)-Llcos(0 e .) 1 e; 
= (_1)]^CG0I+[V(G-)] ]^e^(i^ w /4) IJexp(*^ t ,./4)(-l)"^CGDI i -l^(G0l 

V V* 

=i x(s) ( _ 1} w ; 

with N the number of t> G V(G) such that $ v /2it is congruent to 3 modulo 
4 (and similarly for TV*). By the discrete Gauss-Bonnet formula (recall 
Subsection 12. 3p . 



2ttx(£) = ^(2tt - 0„) + ^(2vr - 

V V* 

This means exactly that N + Af* and ^x(^) have the same parity, so A is 
equal to 1 and the theorem is proved. □ 

4.3. Kac-Ward matrices versus discrete Laplacians. Quite surpris- 
ingly, if the surface £ has genus zero or one, then our Kac-Ward determi- 
nants with critical weights turn out to be proportional to the determinants 
of discrete critical Laplacians. Let us briefly recall the definition and main 
properties of these objects before stating the precise result. 

As first observed by Eckmann [T7j (in a much more general context), the 
Laplace operator on the space of complex valued smooth functions on a Rie- 
mann surface admits a beautifully simple discretization. It is the operator 
A on C°(G; C) = C V W given by 

(Af)(v)= £ x e (f(v)-f(w)), 

e=(v,w) 

for any / G <C V ^ G \ where the sum is over all oriented edges e of the form 
(v, w). Its codimension-one minors are very useful, as they count the number 
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of (weighted) spanning trees in G: this is Kirchhoff 's celebrated matrix tree 
theorem [30]. On the other hand, the determinant of A always vanishes. 

This construction admits a straightforward generalization: 

Definition 4. Let (G,x) be any weighted graph, and (p: tt\(G) — > S 1 any 
representation. The associated discrete Laplacian is the operator A v = 
A*(G, x) on C y ( G ) defined by 

(^f)(v)= xe(f(v)-f(w)<p(e)), 

e=(v,w) 

for any / G C V ( G \ where the sum is over all oriented edges e = (v,w). 

The determinant of this discrete Laplacian has a nice combinatorial in- 
terpretation, which goes back at least to Forman [21J. We include a proof 
here for the sake of completeness. 

Given a graph G, let 3~i (G) denote the set of subgraphs F C G such that 
F spans all the vertices of G, and each connected component T of F has a 
unique cycle (i.e: \V(T)\ = \E(T)\). 

Proposition 4.5 (Forman [21]). For any weighted graph (G,x) and any 
if: tt\(G) -> S 1 , 

detA^(G,*) = Yl nt 2 -^)"^)^)^)' 

Fe3"i(G) TCF 

where the product is over all connected components T of F, x(F) = Y\ e <=F x e> 
and (f(T) = Y\ eG Q T ^(e) with Ct the unique cycle in T endowed with an 
arbitrary orientation. 

Proof. Evaluating directly the determinant of A^(G, x) as a sum over per- 
mutations of V = V(G), we get 

detA^(G,x) = Y (-l) sgnW n A ^)- 

creS(V) vev 

Each permutation a G S(V) decomposes into disjoint cycles, inducing a 
partition V = |_]jVj(°")- A permutation will have a non-zero contribution 
to the determinant only if these orbits correspond to the disjoint union of 
oriented simple closed curves 7 = \_\2Li 7j m ^ A(G) denotes the set of 
such union of curves, we get 

detA*(G,x) = (-1) |7| ^(7^(7)IIE X - 

7eA(G) d^7 e3v 

with x{i) = l\e^r x e and 99(7) = Ueej^)- 

Given any subset W C V, let X(W) denote the set of maps X assigning 
to each v G W an oriented edge X(v) G E with origin v, and set x(X) = 
Ylvew x W(v)- F° r an y "vector field" X G "X(V), let A(X) denote the set 
of disjoint unions of "trajectories" of X: these are oriented simple closed 
curves C in G such that any oriented edge e G C with origin o(e) satisfies 
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X(o(e)) = e. We have 

detA*(G,x)= (-l) l7l v(7)x(7) E x ^ 

7eA(G) xex(v\(vn 7 )) 

= E E (-1) |7| V(7M*) 
xex(V) 7 eApf) 

= e n^-^))^)' 

X£X(V) c 

the product being over all trajectories C of X. (This is the original formula 
of Forman |21j.) The map assigning to each X £ X(V) the subgraph F = 
{X(v)} v< =v is surjective onto the set 3"i(G) defined above, and each F G 
3~i(G) has 2l F l preimages given by the possible choices of orientations of the 
cycle in each of the \F\ connected components of F. The proposition now 
follows from the equality (1 - <p(C))(l - <p(-C)) = 2- tp(C) - ^{C)- 1 . □ 

Let us mention one more fact: if G is a planar isoradial graph, then 
the corresponding critical weights for the discrete Laplacian are given by 
c e = tan(# e ), where 9 e denotes the half-rhombus angle [27j . Therefore, we 
shall refer to A ¥ '(G, c) as the critical discrete Laplacian on G. Finally, we 
shall make the usual abuse of notation and denote by the same letter ip a 
homomorphism vri(S) — > S 1 and the induced homomorphism on vri(G) for 
GcS. 

Theorem 4.6. Let G be a graph isoradially embedded in a flat surface X, 
and let us assume that all cone angles of singularities v 6 V(G) are odd 
multiples of 2ir. Let v e = tan(# e /2) denote the critical weight system on 
GcS, and set c e = tan(# e ). // the genus of T, is zero or one, then for any 

T*(G,v) = (-lf2-xW n iwn detAy(G ' c) ' 

where N is the number of vertices v € V(G) such that "d v /27r is congruent to 
3 modulo 4 and x{G) = |F(G)| — |-E(G)|. On the other hand, the functions 
t^{G,v} and det (G , c) are never proportional if the genus ofT, is greater 
or equal to two. 

Proof. If (p is trivial, then both sides of the equality vanish by Proposi- 
tions HJ] and H2J In particular, the equality holds in the genus zero case, so 
it can be assumed that the genus g of £ is positive. Let F be an element of 
the set 3~(G), that is, a spanning subgraph of G such that no connected com- 
ponent of F is a tree. For any connected component T of F, let Cj, denote 
the corresponding coefficient cl^ = n 7 c9A r (F)(^ ~~ vil))- First note that c!£ 
only depends on the homotopy type of T in E. (Basically, one can deform T 
continuously in £ without changing c!£.) Therefore, it can be assumed that 
T is a wedge of n circles, with n = 1 — x(^) = 1 + 1-^(^)1 — 1^(^)1 > 1- If n 
is greater then 2g, then these n cycles are linearly dependant in ffi(S;Z). 
Via a homotopy of (X,T), one can therefore assume that one of these cy- 
cles is null-homologous, leading to Cj. = 0. If n is equal to 2g, then either 
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these cycles are linearly dependant in Hi(E; Z) and c% vanishes as above, or 
these cycles are independent in homology. In this case, T induces a cellular 
decomposition of E, so dN(T) is the boundary of / faces with 

2 - 2g = X (E) = \V(T)\ - \E(T)\ +f = l-2g + f. 

Therefore, / is equal to 1, so dN(T) is connected, hence null-homologous, 
and vanishes in this case as well. We have proved that c^ F vanishes unless 
each connected component T of F satisfies < — |I/(T)| < 2g — 2. 

In the case of genus 1, this shows that c F vanishes unless F belongs to 
the set 3~i(G). For such an element F, the contribution of each connected 
component T can be easily computed: dN(T) consists of two connected 
components homologous to T and — T, so 

4 = (i- <p(T))0- - = 2 - <p(T) - vp(T)- 1 . 

By Propositions 14.11 and 14. 5^ we now have 

r*(G,v) = C Y, I! (1-^(7))^) 

F&{G) ~/CdN(F) 
FeJi(G) TcF 

= C detA v (G,/i) 

= Cfl y ldet A*(G,c), 

since [x e = ic e . The equality now easily follows from the explicit value for 
the constant C given in Proposition 14.11 

Let us finally assume that the genus of S is greater or equal to two. 
By the argument above, it is enough to show that there is some element 
F G 3~(G) \ 3"i(G) with non-zero coefficient cf p . In particular, we just need 
to find a spanning connected subgraph T C G with |-E(T)| — |V(T)| = 1 and 
such that none of the boundary components of N(T) is trivial in Z). 
Such a T is obtained as follows: choose a spanning tree To C G and add 
two edges of G so that the resulting two cycles are linearly independent in 
i?i(S;Z) but have zero intersection number. (This is possible since E has 
genus g > 1 and G induces a cellular decomposition of E.) The resulting 
graph T satisfies the conditions listed above, and the proof is completed. □ 

Remark 4.7. Boutillier and de Tiliere obtained a similar result in [6, 7J 
(see also [H]): given any graph G isoradially embedded in the flat torus, 
they relate det A^(G, c) with the determinant of some (^-twisted Kasteleyn 
matrix of (Tfj, u), where Tq is the graph associated to G via some variation 
of the Fisher correspondence [19J. Actually, using the methods developed 
in [51 Subsection 4.3], one can show that Corollary 12 of [7] is equivalent 
to our Theorem 14.61 in the flat toric case (that is, when g = 1 and S is 
empty). In our opinion, our approach has several advantages. It is more 
general, as we allow singularities and understand the higher genus while 
Boutillier-de Tiliere only deal with the flat toric case; it is more natural, as 
we work on the same graph G throughout without using one of many possible 
auxiliary graphs Tq', and it is simpler, as our whole proof relies solely on 
Propositions 14. l l and 14.51 The demonstration of Boutillier-de Tiliere, on the 
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other hand, is quite substantial and relies on highly non-trivial results of 
Kenyon [27J and Kenyon-Okounkov |28j . 

We conclude this section with one last remark. Let us assume that G 
is isoradially embedded in the flat torus E = C/A, with A some lattice in 
C, and let G be the corresponding A-periodic planar graph. For n > 1, 
set G n = G/nA. The free energy per fundamental domain of the critical 
Z-invariant Ising model on G is defined by 

f = - lim \\ogZ J (G n ), 

n— »oo n 

with J e = \ log ^ 1 TOsef e ) as explained in Subsection 12.21 The equality 

2 

Z J (G n ) = ( \{ cosh(J e )) n 2 n2 \ v ^\Z(G n ,v) 

e&E{G) 

together with Theorems 13.31 and 14.61 imply 

f = -\V(G)\ l ^P- - ilogdetiA(G,c), 

where log deti A(G, c) = linij^oo -y logdet A A (G n , c) and A is one of the 
three non-trivial spin structures on the torus. (On the torus, spin structures 
are canonically identified with H 1 (Y,; 1,2) , the trivial spin structure has Arf 
invariant 1 and the three others have Arf invariant 0. The trivial spin struc- 
ture does not contribute to the partition function by Proposition 14.11 and 
it is a fact that the three other spin structures will have the same contri- 
bution to the free energy of the model.) The free energy f 1 was computed 
by Baxter in pt], and the normalized determinant by Kenyon in [27] - even 
though the existence of the limit was not proved there. One can check that 
these two results are related in the way displayed above. This gives a reality 
check to our computations, and allows to obtain any of these two results as 
a corollary of the other one. 

4.4. The Kac-Ward determinants as discrete 3-torsions. In this last 
subsection, we wish to relate (in an informal way) the critical Kac-Ward 
determinants with the 3-torsions of the underlying Riemann surface. Let us 
start by briefly recalling the definition of this invariant, in the special case 
relevant to us. 

To any closed N-dimensional complex manifold £ endowed with a unitary 
representation <p>: vri(S) — > U(n) and a Hermitian metric, Ray and Singer 
[35] associate a sequence of numbers T p (£, ip) with p = 0, 1, . . . , N. In the 
case of a Riemann surface (N = 1), the numbers To(T,,p>) and Ti(E,(p) 
coincide, leading to a single invariant T(S, ip). If n = 1, then (p is a unitary 
character (that is, an element of -ff 1 (S; S* 1 )), and it induces a complex line 
bundle L(p) over X. Since S is endowed with a Hermitian metric, one can 
consider the associated Laplacian A on the space of smooth sections of L(p>). 
The d-torsion T(T,,ip) is defined as the square root of the zeta-regularized 
determinant of this Laplacian, that is, 

T(S,^) = exp(-iC'(0)), 
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where £(s) is the zeta function of the Laplacian A. Of course, this depends 
on the choice of Hermitian metric. However, for any two non-trivial char- 
acters ip,f', the ratio T(£, </?)/T(£, ip') is independent of this choice [35j 
Theorem 2.1]. 

Let us come back to the discrete setting. As before, let G be a graph 
isoradially embedded in a flat surface £ with all cone angles at singularities 
in V(G) being odd multiples of 2ir. Recall that the flat metric defines a 
conformal structure on the underlying surface, so that £ is now a Riemann 
surface. Consider the following statement. 

As a function of if £ ff 1 (S;S' 1 ), the Kac-Ward determinant r^(G, v) 
behaves - up to a multiplicative constant - as a discrete version ofT(T,,ip) 2 , 
the square of the corresponding d-torsion. 

We shall not give any proof of this vague statement, not even formulate a 
precise conjecture (although this is very tempting). Instead, we shall simply 
give a list of evidences towards such a statement. 

First of all, recall that given any two non-trivial characters f,f', the ra- 
tio T(£, f)/T(T,, f') only depends on the Riemann surface £ and on f,f', 
but not on the Hermitian metric. Therefore, one expects the correspond- 
ing ratios r^(G, z/)/r </3 (G, u) to exhibit some independence of the choice of 
the graph G embedded in the surface £. As mentioned in the introduc- 
tion, the papers ]12\ 113] provide numerical evidences towards such a claim 
at the scaling limit (in the special case where ip and if' are discrete spin 
structures). Now, our duality result can be understood as a further step 
in that direction. Indeed, it is a direct consequence of Theorem 14.41 that 
the ratios t v (G, v)/^' (G,v) and r^(G*, u*)/t^(G*, u*) actually coincide, 
without taking any scaling limit, for any two non-trivial characters f,f'. 

In the remainder of the discussion, we shall distinguish between the genus 
zero, the genus one, and the higher genus cases. 

The genus zero case. Whenever the character if is trivial, t^{G } v) van- 
ishes by Proposition 14.11 (actually. Lemma 14.21 is enough). Since the kernel 
of the continuous Laplacian has dimension one, T(£, ip) is also equal to zero. 
In particular, t v {G,v) and T(Y,,<p) 2 trivially coincide in the case of genus 
zero. 

The genus one case. In the toric case, Theorem 14.61 shows that r^iG^v) 
is proportional to the most 'natural' discretization of T(£,(/?) 2 , that is, the 
determinant of the critical discrete Laplacian. 

The case of genus g > 1. Let us now consider the case of a Riemann 
surface of genus g > 1. It can be written £ = H/T, where H is the Poincare 
upper half plane and V a discrete subgroup of PSL(2, M). Since E is compact, 
each element 7 E 7q(£) = T acts on H via 7(2) = exp(2p 7 )|5f fl f° r some real 
fixed points Zq, z\ and p 7 > 0. Given any unitary character <p of vri(E) = T, 
the corresponding 3-torsion satisfies 



T(Z,ip) 2 = C-Z r (l,f), 
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with C a constant depending on the Hermitian metric and Z-p (cr, ip) the 
Selberg zeta function defined by 

Zt{?, ¥>) = II II^ 1 - Vb) exp(-2/9 7 (a + A;)))- 

{7} k>0 

Here, the product is over all conjugacy classes of primitive elements of T, that 
is, elements that are not powers in T [35\ Theorem 4.6]. As explained in [36], 
there is a more geometric way to understand this zeta function: conjugacy 
classes of primitive elements of T correspond to primitive closed geodesies 
on S, and 2/3 7 is the length of the corresponding geodesic. Summing up 
(and neglecting the problems of convergence), we have the equality 

T(^f = C [] JJCl-^expC-A:^)), (*) 

7ey(E) k>l 

where denotes the set of primitive closed geodesies in S. 

Let us make the corresponding study on the discrete side. Given an 
oriented closed path 7 on a graph G, we will say that 7 is reduced if it never 
backtracks, that is, if no oriented edge e is immediately followed by the 
oriented edge e. The oriented closed path 7 will be called prime if, when 
viewed as a cyclic word, it cannot be expressed as the product 5 r of a given 
closed path 5 for any r > 2. Finally, we shall denote by 7(G) the set of 
prime reduced oriented closed paths in a graph G. Let us now assume that 
G is isoradially embedded in a flat surface S. As explained in the proof of 
Theorem 13.31 Bass' Theorem |3J implies the equality 

t*(g,u)= n (1-^7)6^(^)1/(7)), M 

where 2/(7) = v e and a 7 is the sum of the angles a(e, e') along 7 (recall 
Figure 

The expression (**) can be understood as a discrete version of (*) in 
the following sense. Primitive closed geodesies on £ are replaced by ho- 
motopically non-trivial reduced prime closed paths in G C S. The factor 
exp (h a j) measures to which extend 7 differs from a straight line (i.e. a ge- 
odesic on E with respect to the flat metric). And finally, the critical weight 
v e = tan(# e /2) plays the role of exp(— £ e ), with l e the length of the edge e. 
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